Interstellar scintillation as the origin of rapid radio vari- 
ability IN THE QUASAR J1819-f3845. 
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Quasars shine brightly due to the liberation of gravitational energy as matter 
falls onto a supermassive black hole in the centre of a galaxy[]T|. Variations 
in the radiation received from active galactic nuclei (AGN) are studied at all 
wavelengths, revealing the tiny dimensions of the region and the processes of 
^ ■ fuelling the black hole. Some AGN are variable at optical and shorter wavelengths, 

O ' 3nd display radio outbursts over years and decades^]. These AGN often also 

show faster variations at radio wavelengths (intraday variability, IDV)[^, Q which 
I have been the subject of much debate[||. The simplest explanation, supported 

O I by a correlation in some sources between the optical (intrinsic) and faster radio 

variations[|, 0], is that the rapid radio variations are intrinsic. However, this 
explanation implies physically difficult brightness temperatures, suggesting that 
2 ! the variations may be due to scattering of the incident radiation in the interstellar 

^ I medium of our Galaxy[|], |^. Here we present results which show unambiguously 

?^ ' that the variations in one extreme case[|T^ are due to interstellar scintillation. 

We also measure the transverse velocity of the scattering material, revealing a 
^ ■ surprising high velocity plasma close to the Solar System. 

H ■ 

\ The most extremely variable extragalactic source known at radio wave- 

lengths is the quasar J1819+3845. This quasar shows variations of factors 
of 4 or more on a timescale of hours, easiest interpreted as scintillation due 
to scattering in the interstellar medium (ISM)||10|. One can think of this 
scattering medium as focusing and defocusing the waves from the source. 



producing a pattern of dark and bright patches[|TT|. As we move through 
these patches, or scintles, we observe a temporal variation in intensity. In 
order to display this effect, the source must be small: J1819-I-3845 is only 
tens of microarcseconds (a few light-months at the source redshift 0.54)|jl0| . 



As the Earth moves through the scintillation pattern, two telescopes At 
different locations should see a difference in the observed time of the intensity 
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Earth viewed from J1 81 9+3845 
12 Jan 




Figure 1: Sketch illustrating the scintles projected onto the solar system. 
The inset shows the Earth and its projected velocity w.r.t. LSR (blue ar- 
row) as seen from J1819+3845. Filled and open symbols show the telescopes 
at start and end of the observation run respectively. As the Earth rotates, 
the value of the 'delay' in arrival of intensity extrema will change during 
the observation due to the changing orientation of the baseline between the 
telescopes with respect to the vector velocity of the scintillation pattern. 
The delay is also dependent on the time of year (due to the Earth's move- 
ment around the Sun) and the relative motion of the Solar System and the 
scattering plasma. 



maxima and minima (fig.|l|). A similar technique was used by Rickett & 
Lang[|l^ who used the decorrelation of signals between telescopes to put 
limits on the size and apparent velocity of the scintles from a pulsar and by 



Jauncey et al[13] for the other known intra-hour variable PKS 0405-387|14 



Compared to other IDV sources, or even PKS 0405-387, the rapidity and 
depth of the variations in J1819-I-3845 make it the perfect candidate for 
measurement of this delay. 

The observations were conducted simultaneously with the NRAO Very 
Large Array (VLA) in New Mexico, and the Westerbork Synthesis Radio 
Telescope (WSRT) in the Netherlands. The source showed large, rapid vari- 
ations on both days (fig. |^). We see that there is a difference in the time at 
which the earliest peaks reach the telescopes (WSRT leads VLA). We also 
see that at the end of the observations this delay has reversed, such that the 
VLA leads WSRT. This delay, and reversal, is repeated on both days. If the 
scattering plasma is stationary with respect to the Local Standard of Rest, 
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the predicted time lag from the motion of the Earth and Sun alone is ~ 
200 sees, with the signal arriving at the VLA first throughout the run. This 
clearly cannot account for the observations, so we used the observations to 
solve for the velocity of the scattering plasma w.r.t. the LSR, transverse to 
the line of sight. 

From a minimization of the total squared flux density difference we cal- 
culate the projected plasma velocity. We obtain for 7 Jan 2001: v^a = 
-32.8 ± 0.5, VDec = 14.3 ± 3 km/s and for 12 Jan 2001: vra = -32.2 ± 0.5, 
VDec = 16.1 ± 1 km/s. It can be seen that we find very good agreement 
between the results obtained from observations of both days. This shows 
that to first order the assumption of a temporally 'frozen' scintillation pat- 
tern is justified. Combining these we calculate the plasma velocity to be 
YjiA = —32.5 ± 0.5, VDec = 15.5 ± 1 km/s (i.e. the plasma moves West and 
North). We also fitted a single time delay to ~ 15 minutes of data around 
the largest maxima, normalising the two lightcurves over the time of interest 
to the same mean flux density. Solutions which pass through these values are 
consistent with the results quoted above. We plot both these results and the 
delays fitted as described above in fig. |^. We can constrain the plasma veloc- 
ity so well because at these values the predicted delay changes very rapidly 
with changes in plasma velocity. The existence of many maxima and minima 
in our data is important as it removes the possibility of the delay, and the 
calculated plasma velocity, being an artefact of flux calibration errors. 

We have previously suggested |T5[] a velocity of the plasma screen in order 
to explain the change in timescale of the variations over a year, but this was 
not a measurement as here, but an interpretation of the data, and the velocity 
was very poorly constrained. For comparison the preliminary result presented 
in ||15[, VRA='i km/s, V£)ec=25 km/s predicts a very different delay: from 
around +60 sec at the beginning of the observations to around +150 sec at the 
end. We will present a combination of both these analyses in a forthcoming 
paper. 

We point out that we are measuring the transverse velocity of the local 
plasma: in general, radial velocities are measured (via spectroscopic tech- 
niques). Also note that whereas scintillation observations of pulsars are sig- 
nificantly affected by their own transverse velocities (up to several hundred 
km/s), the effect of a transverse velocity of (jet features in) AGNs even as 
high as the speed of light is reduced to a fraction of a km/ sec by the distances 
involved. 
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Figure 2: Results of VLA and WSRT observations from 7 and 12 January 
2001. Data were averaged over 30 sec. The time difference in tlie arrival of 
the signals can be seen as a horizontal separation of the blue and red lines. 
The vertical axis of the difference signal before and after fitting has been 
expanded by a factor of 2.5. A total of 11 VLA antenna were used for the 
observations described here. The central frequency at both telescopes was 
4875MHz, and the bandwidths similar (80 and 100 MHz). VLA observa- 
tions of a nearby calibrator source (J1821+3945) are shown in the top panel 
of each plot. Slightly different observing strategies were used on each run. 
On 7 January 2001 we observed the source continuously, with a sub-set of 
antennae observing the calibrator. On 12 January 2001 we used the same 
antennae to nod between the source and the calibrator, observing the latter 
for 3 minutes every ~ 35 minutes. On both dates the observations described 
above were bracketed with observations of the calibrator and a flux density 
calibrator (3C286). Elevation (and other time dependent) gain variations 
at the VLA were determined using the calibrator J1821-I-3945 were found to 
be up to 3%. The data was corrected for these effects (see below). Obser- 
vations in good weather at WSRT typically yield 1% flux density reliability. 
This is applicable to all but the first hour of the 7 January 2001 run, when 
light rain caused an estimated extra 1-2% extinction. The error we quote 
on the result is obtained using a nuiA)er of different determinations of the 
absolute flux density scale, and acceptable changes throughout the run, in- 
cluding: no temporal changes in gain corrections, a linear change, and gain 
corrections determined by observations of the calibrator, which is assumed 
stable. Thermal noise in 30 sec is ImJy (WSRT) and 0.6 mJy (VLA). 
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Figure 3: The calculated delays as a function of observing time for both days. The 
delays are calculated as a function of time from the projected baseline and velocity 
vector. We use a circular Earth orbit and assume that the Sun is moving at 19.7 
km/s towards 18h, 30° (J2000). As the direction of the Sun's motion happens 
to be only 10° from J1819+3845, errors caused by its uncertainty are negligible. 
For a wide range of plasma velocities we calculated the total squared flux density 
difference, after time shifting the data according to the calculated delay. Mini- 
mizing this quantity gives the plasma velocity. The associated delays are plotted 
as the diagonal lines. The results from a given calibration scheme produce errors 
of less than 1 km/sec, but the final error is dominated by calibration uncertainty. 
The error we quote is that obtained using a number of different determinations 
of the absolute flux density scale, and acceptable changes throughout the run. In 
addition, independent fits to portions of the data are shown with 3a error bars. 
The data from times with large delays are shown in the insets (red and dark blue 
as before, yellow showing the WSRT data after fitting). 
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The determination of a peculiar velocity argues that the scattering struc- 
ture is discrete. This measurement of the transverse velocity of the ionised 
scattering material cannot be directly compared to radial velocity measure- 
ments of HI, but we note that there is no anomalous feature known in this 
direction |]16|. A nearby pulsar J1813-I-4013 has an unremarkable dispersion 



and rotation measures (42 cm pc and 47radm )|r?], although the medium 



in front of J1819+3845 may not extend to this line of sight. We have pre- 
viously estimated [jlO[ that this scattering region is located at around 20 pc 



from the Sun. The scatterer may be related to the Local Bubble||18[, and if 



so, for the reasons given above, is more likely to be a discrete structure at 
the boundary than due to the putative intervening hot gas. 

A small difference in the flux densities observed at the two telescopes 
(after delay correction) as the telescopes pass through different regions of 
the scintles is expected. The size of the scintles in the direction of motion 
is ~ 4000 sec x 25 km/sec = lO^km. If the scintles were of similar size per- 
pendicular to this, we can estimate the expected difference in the intensities 
on a typical baseline of 6000 km by considering the flux density differences 
between the two curves if they are shifted 6000/10^ x 4000 = 240sec. This 
gives a potential 5-10% flux density difference (after delay correction) be- 
tween the two curves. However, differences are less than 2% on both days 
(fig.^, suggesting the scintillation pattern is highly anisotropic, in agreement 



with preliminary analysis of the variation of the timescale over a year[^ . 

To account for these observations there must be an unusually high fraction 
of turbulent electrons, if the electron density is that of the typical extended 
ISM. In comparison to the conditions found in highly turbulent, ionised re- 



gions, e.g. around OB stars p9|, pO| , |21| however, the required conditions are 
very modest. Interestingly, the AO star a Lyrae (Vega), which at Spc is 
within the estimated distance range for the screen, is at a projected angular 
separation of J1819+3845 of 3°, which is only 0.4 pc. Whether the origin of 
the turbulent plasma could be related to mass loss from this star, or a feature 
of the so-called Local Fluff p2[, remains to be investigated. 

Our observations prove beyond doubt that the variations in J1819-I-3845, 
the most extreme example of an IDV source, are due to scintillation. Yet, as 
this source must be extremely small, it could easily be intrinsically variable on 
time scales of a few months (using simple light-travel time arguments). Using 
relativistic outflow arguments its is not difficult to bring this down to weeks. 
Intrinsic variations on these time scales would be very hard to separate from 
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scintillation induced variations. We have evidence however that the source 
is remarkably stable over a period of many months with only a slow increase 
in brightness over the last two years (Dennett-Thorpe & de Bruyn, in prep). 
Scintillation in known variable sources and those with slower variations will 
be more difficult to prove directly. However, the few percent variations seen 
on day scales in many IDV|2^ is likely to be due to scattering of structures 
of tens of microarcseconds through more typical ISMp4[. This implies that 
most AGN have brightness temperatures no higher than 10^^ — lO^^K, not 
10^^ — 10^°K, as would otherwise be required. Perhaps the most serious 
argument against an intrinsic explanation of the IDV phenomenon, which 
requires extremely small angular sizes, now becomes how to avoid scintillation 
induced variations dominating all the observed variability. 

The lack of correlation of strength or timescale of variations with fre- 
quency or galactic latitude is sometimes cited as evidence against ISS as 
cause for IDV [p5|] . However as the variations are affected by the intrinsic size 
of the source, the location and motion of the effective scatterer as well as 
the season of the observations, such correlations may easily be washed out. 
Although the annual modulation in the variations of B 0917+62 have been 
explained in terms of scintillation through the general Galactic ISMp6|, |27| , 
J1819+3845 demonstrates that one cannot assume that the scattering ef- 
fectively occurs throughout this medium. It is important to note that the 
scattering of resolved AGN will be weighted by the material closer to the 
observer. 

The delay in arrival times of the intensity variations at two widely sepa- 
rated radio telescopes is a neat proof that a propagation effect causes the ex- 
treme variations in this one quasar. The small size of the quasar J1819+3845 
is probably related to the unusual rising radio spectrum. The absence of 
many sources which scintillate in this manner is presumably due to a com- 
bination of an unusually small, stable source, and an atypically strong, close 
scatterer. We nonetheless expect that more similarly extreme IDVs will be 
discovered. 
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